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We identify stationary distributions of generalized Fleming- Viot pro- 
cesses with jump mechanisms specified by certain beta laws together 

■ with a parameter measure. Each of these distributions is obtained from 
Oh . normalized stable random measures after a suitable biased transforma- 
r S^ tion followed by mixing by the law of a Dirichlet random measure with 

the same parameter measure. The calculations are based primarily on 
the well-known relationship to measure- valued branching processes with 
l— 'I immigration. 

(N : 

^ ; 1 Introduction 

In the study of population genetics models, it is of great importance to identify their 
stationary distributions. Such identifications provide us with basic information of 
possible equilibria of the models and are needed prior to quantitative discussions on 
statistical inference. Since [5], [14] and [1], theory of generalized Fleming- Viot pro- 
cesses has served as a new area to be cultivated and has been developed considerably. 

■ (See [3] for an exposition.) In view of such progress, it seems that we are in a position 
to explore the aforementioned problems for some appropriate subclass of those mod- 
els. In this respect, it would be natural to think of the one-dimensional Wright-Fisher 
diffusion with mutation as a prototype. This celebrated process is prescribed by its 
generator 

A:=^x(l-x)-^ + ^[ Cl (l-x)-c 2 x}-^, are [0,1], (1.1) 

where c\ and c 2 are positive constants interpreted as mutation rates. The stationary 
distribution is a beta distribution 

B CuC2 (dx) : = £^±^1^-1(1 _ %r -Hx, (1.2) 

r ci r c 2 
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where T(-) is the gamma function. In addition, the process associated with (1.1) 
admits an infinite-dimensional generalization known as the Fleming- Viot process with 
parent-independent mutation, whose stationary distribution is identified with the law 
of a Dirichlet random measure. 

In the present paper we consider a problem of finding a class of generalized 
Fleming- Viot processes whose stationary distributions can be identified. As far as 
the first term on the right side of (1.1) is concerned, its jump- type version has been 
discussed in population genetics as the generator of a model with 'occasional extreme 
reproduction'. (See Section 1.2 of [3] for a comprehensive account.) We addition- 
ally need to look for an appropriate modification of the second term, which should 
correspond to a generalization of the mutation mechanism. With these situations in 
mind, our problems can be described as follows. 

(I) By modifying both mechanisms of reproduction and mutation, find a jump pro- 
cess on [0,1] whose generator extends (1.1) and whose stationary distribution can be 
identified. 

(II) Establish an analogous generalization for the Fleming- Viot process with parent- 
independent mutation. 

Since these problems are rather vague, it may be worth showing now the generator 
we will believe to give an 'answer' to (I). For each a G (0, 1), define an operator A a 



where G are smooth functions on [0, 1]. Observe that A a G(x) — > AG{x) as a f 1. 
It should be noted that A a is a one-dimensional version of the generator of the 
process studied in [2] if c\ = c 2 = 0. See also [12] and [13]. The reader, however, 
is cautioned that our notation a is in conflict with that of these papers, in which 
a plays the same role as a + 1 in our notation. (We adopt such notation in order 
for the formulae below to be simpler.) The constant C\ (resp. c 2 ) in (1.3) can be 
interpreted as the rate of 'simultaneous mutation' from one type to the other type 
and a proportion u of the individuals with that type, which are supposed to have 
the frequency 1 — x (resp. x) in the population, are involved in this 'mutation' event 
with intensity B\^ aM {du) / ((a + l)u). (Note that (1 — u)x + u = x + u{l — x).) As 
will be seen in Proposition 3.1 below for more general case, the closure of (1.3) with 
a suitable domain generates a Feller semigroup on C([0, 1]), and our main concern 
is the equilibrium state of the associated Markov process. It will be shown in the 
forthcoming section that a unique stationary distribution of the process governed by 
(1.3) is identified with 



by 



A a G(x) 




[xG((l - u)x + u) + (1 - x)G{{\ - u)x) - G(x)} 




r 1 B^g^du) 
o (a + l)u 



[ciG((l - u)x + u) + c 2 G({l - u)x) - (ci + c 2 )G(x)\ (1.3) 




Y 1 +Y 2 



Y 1 



G dx 



(1.4) 
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where E ajV denotes the expectation with respect to (Yi,Y 2 ) with law determined by 
\ogE a: y{e- XlYl - X2Y *] = -j/A? - (1 - y)\% (Ai, A 2 > 0). Again we see that (1.4) with 
a — 1 reduces to (1.2). 

One might think that (1.3) appears very special among possible generalizations 
of (1.1). In fact, such restriction has resulted from our strategy, which is described 
in the following. It is well-known [20] that the Fleming- Viot process with parent- 
independent mutation can be obtained by way of a normalization and a random time 
change from a measure-valued branching diffusion with immigration. (See also [10] 
and [18] for related works.) An extension of this significant result was shown in [2] for 
a class of generalized Fleming- Viot processes, one-dimensional version of which corre- 
sponds to (1.3) with ci = c 2 = 0. Moreover, [2] proved that such a jump mechanism 
is necessary for a generalized Fleming- Viot process to have the above mentioned link 
to a measure- valued branching process with immigration (henceforth MBI-process). 
Recently, [13] showed essentially that the second term of (1.3) is required when we 
additionally take a generalization of the mutation mechanism into account. Our ar- 
gument will be crucially based on this kind of relationship between the generalized 
Fleming- Viot process associated with a natural generalization of (1.3) and a certain 
ergodic MBI-process. That relationship can be reformulated as a factorization result 
on the level of generators and hence is expected to yield also an explicit connection 
between stationary distributions. In principle, the problems (I) and (II) can be con- 
sidered in a unified way. Nevertheless, we shall discuss (I) and (II) separately. This is 
mainly because the factorization identity will turn out to yield a correct answer only 
for certain restricted cases and in one-dimension one can avoid to use it by taking 
an analytic approach instead (although this does not reveal clearly the mathematical 
structure underlying). 

The organization of this paper is as follows. Section 2 is devoted to derivation of 
(1.4) by purely analytic argument. Exploiting the relationship to MBI-processes, we 
show in Section 3 that the above mentioned answer to (I) has a natural generalization 
which settles (II). The irreversibility of the processes we consider is discussed in 
Section 4. 

2 The one-dimensional model 

Let < a < 1, Ci>0 and c 2 > be given. The purpose of this section is to show 
that (1.4) is a unique stationary distribution of the process with generator (1.3). 
Analytically we shall prove that a probability measure P on [0, 1] satisfying 



is uniquely identified with (1.4). Actual starting point of the calculations below is 




for all G(x) = <f n ( x ) '■= x n with n — 1,2,... 



(2.1) 




for all G(x) = G t (x) := (1 + tx) 



1 with t > 0. 



(2.2) 
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The equivalence of (2.1) and (2.2) is a consequence of uniform estimates 

\A a <p n (x)\< (i + ^±^W n = l,2,..., 
V a + 1 / 

which can be shown by observing that 

ci((l - u)x + u) n + c 2 ((l - u)x) n - (ci + c 2 )a;" 
= ci [((1 - u)x + «)" - ((1 - u)x + wa;) n ] + c 2 x n [(1 - u) n - ((1 - u) + u) n ] 



ci £ ( " ] (1 - «) n " V" V(l - :r fe ) - c 2 x" ECJK 1 - 



fc=i 



ei;i(i- 



k=l 



k, 



u) n k u k 



u) n k u k 



k=l 



Cl x n - k - (ci + c 2 )x r ' 



(2.3) 



= «E(JJ(i-«r« w 

and in particular 



cix n fc - (ci + c 2 K' 



x((l - u)z + u) n + (1 - x)((l - w)x) n - x n 



= Yl [ J (1 - «) n "V(a; n - fc+1 - -" 



x n ) 



fc=2 



Indeed, these bounds ensure that the function 



-1 oo x 

^ / A a G t (x)P(dx) = J2(-t) n A aVn (x)P(dx) 
Jo n=l Jo 

is real analytic at least for —1/2 < t < 1/2. We prepare a simple lemma in order to 
calculate A a G t - 



Lemma 2.1 Assume that b > and a + b > 0. 

(%) It /jo/ds £/ia£ /or any #i > and 9 2 > 



Jo 



Be ud2 (du) 



= (a + b)-H- e \ 



{au + 

(^zj In addition, suppose that a' ^ a and a' + b > 0. T/ien 
5i- a ,i+a(dw) 1 



(2.4) 



I 1 
Jo 



(au + b)(a'u + b) a (a - a')b 1+a 



a + b) a - (a' + b) c 



(2.5) 



(2.4) is a one-dimensional version of the formula due to [4], which is sometimes 
referred to as the Markov-Krein identity. (See e.g. [22] or (3.6) below.) We will give 
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a self-contained proof based essentially on the well-known relationship between beta 
and gamma laws. 

Proof of Lemma 2.1. The proof of (2.4) is simply done by noting that 



(a + b)-H- 9 * = ( 
Jo 



dzi ,fli-l r -(a+W*i f°° dz ^ 7 e 2 -l p -bz 2 

o mf 1 J r(e 2 f 2 



and then by change of variables to u := z 1 /(zi + z 2 ),v := zi + z 2 . (2.5) can be deduced 
from (2.4) with 9 1 — 1 — a and 9 2 = a since Bi_ al+a {du) = Bi_ aa {du){l —u)/a and 



1 — u 



a + b a' + b 



(au + b)(a'u + b) (a — a')b \au + b a'u + b / 



We proceed to calculate A a G t . 
Lemma 2.2 For any t > and x e [0, 1], 

^G,W= f . (1 + f) °- 1 .' (1 -'> 

v ; a (l + tx) 2+Q a + 1 

Proof. By straightforward calculations 



* ci(l -x)(l +t) a ~ 1 -c 2 x 



(l + tx) 1+a 



ciG t ((l - w)x + «) + c 2 G t ((l - - (ci + c 2 )G t {x) 



tu 



1 + tx 



Ci(l - x) 



c 2 x 



1 + t(l - u)a; + tit l + t(l-u)x 



Replacing ci and c 2 by x and 1 — x, respectively, we get 

xG t ((l - u)x + u) + (1 - x)G t ((l - u)x) - G t (ar) 
t 2 u 2 x(l - x) 1 



(2.6) 



1 + tx (1 + t(l - u)x + tu)(l + t(l - u)x)' 
Plugging these equalities into (1.3) with G — G t and then applying Lemma 2.1 yield 

t 2 x(l — x) f 1 Bi- a! i +a (du) 



A a G t {x) 



-a,l+a \ 

1 + tx Jo (1 + t(l - u)x + tu)(l + t(l - u)x) 
t fi B^Adu) 

K 'Jo l + t(l-u)x + tu 
" l B L ^ a (du) 



(a + l)(l + tx) 
t 



■ c 2 x 



(a + l)(l + tx) ~* Jo l + t(l-u)x 
t 2 x(l-x) 1 



1 + tx at{l + tx) l+a 



[(1+tr-l] 



(a + l)(l + tx) 



ci(l-x) 



c 2 x 



(1 + t) 1 - Q (l +tx) a (l + tx) a 
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which equals the right side of (2.6). 

Next, we are going to characterize stationary distributions P in terms of 

P{dx) 



*■<«>== jfrr 



t > 0, 



(1 + tx) a) 

which is a variant of the generalized Stieltjes transform of order a. 



(2.7) 



Proposition 2.3 A probability measure P on [0, 1] is a stationary distribution of the 
process associated with (1.3) if and only if S a defined by (2.7) satisfies for all t > 



a 



(2.8) 



+ 



Ci + 1 + -) ((1 + t) a - 1) + ci + c 2 

a J 



S' a (t) + ac 1 (l + t) a - 1 S a (t) = 0. 



Proof. By virtue of Theorem 9.17 in Chapter 4 of [7], P is a stationary distribution 
of the process associated with A a if and only if (2.1) (or (2.2)) holds. By Lemma 2.2 
(2.2) now reads for all t > 



;i + t)°-i r 1 x(i-x 



a 



o (l + tx) 2 +' 



-P{dx) 



This equation becomes (2.8) by substituting the equalities 



1 x(l — x) 
o (1 + tx) 2+a 

1 1 -x 



o (l + tx) 1+a 



P(dx) = — S'Jt) + S„(t) 
a 



and 



1 (1 



(1 + tx) 1+a 
all of which are verified easily. 



P(dx) = —S' a (t), 



a 



We now derive (1.4) as the unique stationary distribution we are looking for. 
Recall that for each y e (0, 1) we denote by E a:V the expectation with respect to the 
two-dimensional random variable (Yi,F 2 ) with joint law determined by 

E a Je- x ^- x ^} = e-^-a-w)*?, Ax, A 2 > 0. 

By using = r(a)" 1 / oo (iOT a ' 1 e _ri (t > 0) and Fubini's theorem, observe that 



E, 



(tY 1 + Y 2 )- a 



dvv"- 1 exp [-y(vt) a - (1 - y)v a ] 
1 



re 

Via)' 1 

Jo 

1 

r(a + 1) ' 1 + (t a - l)y 
6 



(2.9) 



for t > 0. In particular, E a%y \{Y\ + Y 2 ) a ] = 1/T(a + 1) and hence 



(dx) = F(a + l) / B CuC2 {dy)E a , 

J 



(Y 1 + Y 2 y a ; 



Y 1 



Y\ + Y 2 



(2.10) 



defines a probability measure on [0, 1]. Although for each y e (0, 1) an expression of 
the distribution function 



[0, l]3x^T(a + l)E aji 



(Y 1 + Y 2 )- a ; 



Y 1 



Y! + Y 2 



< x 



is given as the formula (3.2) in [23], i.e. 

(1 — y)(x — u) a ~ l u a du 



sin cot r x 

71 Jo (T 



y) 2 u 2a + y 2 {\ - u) 2a + 2y(l - y)u a {\ - u) a cos cm' 



we do not have any explicit form concerning P a ,( Cl ,c 2 ) except the case c\ + c 2 = 1. 
(See Remark (ii) at the end of this section.) 

The main result of this section is the following. 

Theorem 2.4 The process associated with (1.3) has a unique stationary distribution, 
which coincides with P Q) ( Cl ,c 2 )- 

Proof. Notice that the existence of a stationary distribution follows from compactness 
of the state space [0,1]. (See e.g. Remark 9.4 in Chapter 4 of [7].) Let P be 
an arbitrary stationary distribution of the process associated with (1.3) and S a be 
defined by (2.7). Put 

T a (u) = S a ((l + uy/ a -i) 
for u > 0. Setting t = (1 + uf/ a - 1 or u = (1 + t) a - 1, observe that for u > 



TM = -(l + ^- l S' a (t) 
a 



and 



t'M = I(I_i)(i + „)i-X(0 + [-(i + «)H 2 ^(0 

a \a J la J 

= (± - l) (1 + u)-'T' a {u) + 1(1 + u)i- 2 S>>{t). 

Hence S' a (t) = a(l + u) l -^T' a (u) and 

S''(t) = a 2 (l + u) 2 -i 



T^u)-[--l)(l + u)- 1 T^u) 



Also, (2.8) can be rewritten as 



^(1 + + 



ci + H — I w + ci + c 2 



^(t) + aci(l + «) 1 ^ 1 S a (t) = 0. 



From these preliminary observations, it is direct to see that the equation (2.8) is 
transformed into a hypergeometric equation of the form 

u(l + u)T^(u) + [(c 1 + c 2 ) + (c 1 + 2)u]T , a (u) + c 1 T a (u) = 0, u > 0. (2.11) 

Clearly T a (0) = S a (0) = 1. In addition, 

^(0) = S' a (0)/a = - f P(dx)x = - Cl /( Cl + c 2 ), 

where the last equality follows from (2.1) with n — 1. These facts together imply 
that 

4 B CUC2 (dy) 



T a {u) 



or 



o 1 + uy 

B Cl , C2 (dy) 



u > 



/o i + {(i + Q -i}y' 

(See e.g. Sections 7.2 and 9.1 in [16].) Combining this with 
1 , f 1 1 



* > 0. 



r ' a + 1 »i a 



which is immediate from (2.9), we arrive at 



+ t x y Ea > y 



(Y 1 + Y 2 )- a ; 



Y 1 + Y 2 



G dx 



(ci,c 2 ) 



(dx) 



+ tx) a ' 



t > 



(2.12) 



in view of (2.10). Therefore, we conclude that P = P a ,(ci,c 2 ) an d the proof of Theorem 
2.4 is complete. ■ 

Remarks, (i) In the case where c\ + c 2 > 1, an alternative expression for P a ,(a,c 2 ) 
exists: 



P a ,(c u c 2 )(dx) = F(a + l)(ci + c 2 - 1)£ 



(Z 1 + Z 2 )- Q ; 



G dx 



(2.13) 



^1 + ^2 

where Z\ and Z 2 are independent random variables with Laplace transforms 

E[e- xz >] = exp [-a log(l + A Q )] , A > 0. (2.14) 



This reflects the fact that the solution to (2.11) with the same initial conditions 
T a (0) = 1 and T^(0) = —C\j{c.\ + c 2 ) admits another integral expression of the form 



T a {u) 

and accordingly by (2.12) 



1 B 1 



,ci+c 2 — 1 



(dy) 



o (1 + uy) Cl 



u>0 



f l Pg,(cuc 2 )(dx) = r 1 
Jo (l + tx) a Jo 



+C2-1 



[i + {(i + a -i}y] C1 ' 
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t > 0. 



(2.15) 



On the other hand, it is not difficult to show that (2.15) with P a ,( ci ,c 2 ) i n place of 
Pa,(ci,c 2 ) holds, too. In fact, we prove in Lemma 3.5 below a generalization of the 
coincidence (2.13) in the setting of random measures. Also, the role of Z 1 and Z 2 
will be made clear in connection with branching processes with immigration related 
closely to the process generated by (1.3). (Compare (2.14) with (3.9) below.) 

(ii) It will be shown in the Remark after Lemma 3.5 below that -P a ,(ci,c 2 ) = B acit0lC2 
holds whenever c\ + c 2 = 1. At least at a formal level, this would be seen by letting 
Ci + c 2 j. 1 in (2.15) and then by making use of (2.4). 

(iii) In contrast with the case of the Wright-Fisher diffusion mentioned in the Intro- 
duction, P a ,( Cl ,c 2 ) with < a < 1 is not a reversible distribution for the generator 
(1.3) at least in case C\ ^ c 2 . This will be seen in Section 4. 



3 The measure-valued process case 

The main subject of this section is an extension of Theorem 2.4 to a class of gen- 
eralized Fleming- Viot processes. But the strategy will be different from that in the 
previous section, and so an alternative proof of Theorem 2.4 will be given as a by- 
product. To discuss in the setting of measure- valued processes, we need new notation. 
Let E be a compact metric space having at least two distinct points and C(E) (resp. 
B + {E)) the set of continuous (resp. non-negative, bounded Borel) functions on E. 
Define Ai(E) to be the totality of finite Borel measures on E, and we equip M.(E) 
with the weak topology. Denote by M.(E)° the set of non-null elements of Ai(E). 
The set J\A.\{E) of Borel probability measures on E is regarded as a subspace of 
M.(E). We also use the notation (r/, f) to stand for the integral of a function / with 
respect a measure 77. For each r G E, let S r denote the delta distribution at r. Given 
a probability measure Q, we write also E Q [-] for the expectation with respect to Q. 

Let < a < 1 and m G J\4(E) be given. We shall discuss in this section an 
J\4\{E)- valued Markov process associated with 

A a ,m^) (3.1) 

:= f / H(dr) M(l - u)n + u5 r ) - $(„)] 

JO U z JE 

+ f B j- a ^ u) f m(dr) [$((1 - u)im + u5 r ) - $(/,)] , im G M^E), 
Jo (a + l)u Je 

where $ belongs to the class J-'i of functions of the form $/(/i) := (/x® n , /) for some 
positive integer n and / G C(E n ). (3.1) shows clearly that A a ,m satisfies the positive 
maximum principle and hence is dissipative. (See Lemma 2.1 in Chapter 4 of [7].) 
We begin by seeing that A atm defines a Markov process on M.\{E) in an appropriate 
sense. For this purpose, we need an expression for A a , m &f with / G C(E n ). Set 
(a) b = T(a + b)/T(a) for a > and b > 0, and let | • | stand for the cardinality. It 
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holds that for any 9 > and v G M\(E) 



E (1 a)k -f a , + l)n ~ k E «/^ ©5 n V) - S/M) 



r(n) 



fc=2 " V"7 /:|/|=fc 

" (1-aWa) 



+*E 

k=l 



In—k 



{a + l)r(n) 



E «^,SW/) -$/(//)), (3.2) 



I:\I\=k 



where I are nonempty subsets of {1, . . . , n}, 0^ : C(E n ) — > C(E n ) is defined by let- 
ting be the function obtained from / by replacing all the variables with i £ I 
by r m i n j and : C(E n ) — >■ C(E n ) is defined by letting S^/ be the function ob- 
tained from / by replacing all the variables with i G / by r and then by integrating 
with respect to v{dr). (For a degenerate z/ (3.2) is a special case of the corresponding 
expression found in the proof of Lemma 11 in [12].) (3.2) can be deduced from the 
following identities (cf. (2.3)) among signed measures on E n : 



((1 - u)fi(dri) + u5 r (dri)) - fi(dri) 



i=i 



i=i 



((1 — u)n(dri) + u5 r (dri)) — 6?) ((1 — u)ji{dri) + ufj,(dri)) 



i=i 



i=i 



= E0((!-«)M^) 

/^0 jgj 



(uS r (dri)) - (ufi(dri)) 



iei 



iei 



IjLQ jg/ 



i S r (dri 



i //(dr;) 



Lie/ 



ie/ 



As for the Fleming- Viot process with parent-independent mutation, the result corre- 
sponding to the next proposition is a special case of Theorem 3.4 in [8]. 

Proposition 3.1 For each m G M.{E) the closure of A a , m defined on T\ generates 
a Feller semigroup on C{M.\{E)). 

Proof. Let 9 > and v G M.\(E) be such that m = 9v. We simply mimic the 
proof of Theorem 3.4 in [8]. In particular, the Hille-Yosida theorem (Theorem 2.2 
in Chapter 4 of [7]) will be applied. Let n be an arbitrary positive integer. Rewrite 
(3.2) as 

Aafi&ffjx) = (n® n , e (n) f) + e(^ n , s<»>/> - Cn (a, e)$ f (ii), 

where Q^ n \ : C(E n ) — > C(E n ) and c n (a,9) are respectively the non-negative 
operators and the positive constant defined implicitly by the above equation combined 
with (3.2). Let A > be arbitrary. Given g G C(E n ), define 



h = (A + c n (a, 9))- 1 ]T [(A + c n (a, 0))' 1 (e< n > + 9~^) 



k=0 



9- 



10 



Then h G C(E n ) since the operator norm of + 8EW equals c n (a, 0). Moreover, 



so (A — A a fiv)^h — This implies that the range of A — A a ,ev contains T\, which 
is dense in C(M\(E)). The rest of the proof is the same as that of Theorem 3.4 in 



For simplicity, we call the «4 Qim -process the Markov process governed by A a , m i n 
the sense of Proposition 3.1. This process is a natural generalization of the process 
generated by (1.3) in the following sense. Suppose that E consists of two points, say 
T\ and r 2 , set m = C\5 ri + C2<5 r2 , and let {X(t) : t > 0} be the process generated by 
(1.3). Then, verifying the identity A a , m &(n) = A a G(x) for /i = x5 ri + (1 — x)5 r , 2 and 
$(/i) = G(x), we see that the process {X(t)8 ri + (1 — X(t))8 r2 : t > 0} defines an 
^4 ajm -process. We note that [13] discusses the case where E = [0, 1] and m = cS for 
some c > 0. 

We could also establish the well-posedness of the martingale problem for A a , m 
by modifying some existing arguments. More precisely, the existence could be shown 
through a limit theorem for suitably generalized Moran particle systems by modifying 
those considered in the proof of Theorem 2.1 (especially (2.2)) of [14], which took ac- 
count of the jump mechanism describing simultaneous reproduction (sampling) only, 
so that simultaneous movement (mutation) of particles to a random location (type) 
distributed according to m{dr)/m{E) is allowed. The uniqueness would follow by the 
duality argument employing a function- valued process as in the proof of Theorem 2.1 
of [14]. Its possible transitions and the associated transition rates are found in (3.2). 
The duality would be useful in discussing (weak) ergodicity of the ^4 Qjm -process. (See 
e.g. Theorem 5.2 in [8] for such a result in the Fleming- Viot process case.) 

The following argument is based primarily on the relationship between the A a , m - 
process and a suitable MBI-process, which takes values in Ai(E). More precisely, 
the generator, say £ a , m , of the latter will be chosen so that for some constant C > 



where ^(rj) = Q(i](E)~ 1 i]) and <£> is in the linear span Tq of functions of the form 
/i i — y ([/,, /i) • • • (/x, f n ) with fi G C(E), i = 1, • • • , n and n being a positive integer. 
In the case of the Fleming- Viot process (which corresponds to a — 1 formally), 
such a relation is well-known. For instance, it played a key role in [20]. As for the 
generalized Fleming- Viot process, factorizations of the form (3.3) have been shown 
in [2] for m = (the null measure) and in [13] for degenerate measures m. From 
now on, suppose that m G M.(E)° . To exploit (3.3) in the study of stationary 
distributions, we further require the MBI-process associated with £ Qjm to be ergodic, 
i.e., to have a unique stationary distribution, say Q a ,m, supported on A4(E)°. Once 
these requirements are fulfilled, (3.3) suggests that 



(\ + c n (a,e))h-(e^ + ezW)h = g, 



[8]. 





(3.3) 



P Q , m (-) := E Qa ' m \r](E)~ a ; i](Ey 1 i] G •] / E Q ^ \q{E) 



(3.4) 
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would give a stationary distribution of the *4 Qim -process provided that r](E)~ a is 
integrable with respect to Q a ^ m - This conditional answer may be modified to be a 
general one, which must be consistent with the one-dimensional result (1.4). 

To describe the answer, we need both the a-stable random measure with parame- 
ter measure m and the Dirichlet random measure with parameter measure m, whose 
laws on A4(E)° and Ai\(E) are denoted by Q a , m and V m , respectively. These infinite 
dimensional laws are determined uniquely by the identities 

/ Q a , m {dri)e- {vJ) = e~ {m ' n (3.5) 
Jm{e)° 

and 

/ V m {dn)(n, 1 + f)- m{E) = e-Wog(i + /)> ; (36) 

JMi(E) 

where / G B + (E) is arbitrary. A random measure with law Q a , m is constructed from 
a Poisson random measure on (0, oo) x E. (See also Definition 6 in [22].) Observe from 
(3.5) that E Qa - m [r](E)~ a ] = l/(m{E)V{a + 1)). As in [11], V m is defined originally 
to be the law of a random measure whose arbitrary finite-dimensional distributions 
are Dirichlet distributions with parameters specified by m. The useful identity (3.6) 
is due to [4] and reduces to (2.4) in one-dimension. We now state the main result of 
this paper. 

Theorem 3.2 For any m G M.(E)° , the A a , m -process has a unique stationary dis- 
tribution, which is identified with 

P Q>m (.) := r(a + 1) / V m {dn)Eft°>* UE)' a ; ^E)'^ G •] . (3.7) 

To illustrate, consider the trivial case where m = 95 r for some 9 > and r G E. Then 
it is verified easily that P a , m concentrates at S r G Aii(E), and this is consistent with 
the equality A a ,m^(S r ) = in that case. Also, for every m G Ai(E)°, we note that 
P a ,m — > Dm as ct f 1 since by (3.5) Q a ,^ converges weakly to the delta distribution 
at fi for each \i G M.\{E). 

The proof of Theorem 3.2 will be divided into three steps. As mentioned earlier, 
we first find an ergodic MBI-process whose generator satisfies (3.3) and show, under 
necessary integrability condition, that P a ,m i n (3-4) gives a stationary distribution of 
the A^m-process. (In fact, the condition will turn out to be that m(E) > 1. This 
motivates us to make a reparametrization m =: 6u with 9 > and v G M.\{E\.) 
Second, for each v G A4\(E), we prove that P a ,eu = Pa,8u for any 9 > 1. As the last 
step, we extend stationarity of P a fi v with respect to A a fiv to all 9 > by interpreting 
the condition of stationarity as certain recursion equations among moment measures 
which are seen to be real analytic in 9 > 0. Also, the recursion equations will be 
shown to yield uniqueness of the stationary distribution. 

For the first step, we prove in the next proposition that the MBI-process with the 
following generator is the desired one: 
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a + 1 f°° dz r 



5^ 

#(?7 + z5 r ) - #(77) - <^(r) 



1, 

W'! - ' 

ct 077 



4 



where \P is in the class J 7 of functions of the form 77 1->- F((rj, . . . , (77, /„)) for some 
F G C 6 2 (R n ), /j G C(E) and a positive integer n, and g(r) = f (77 + e£ r )| . Up 
to this first order differential term, the operator (3.8) for E = [0, 1] and m = c5 with 
c > is the same as the one discussed in Lemma 17 of [13], in which the factorization 
(3.3) has been proved. Thus, our main observation in the next proposition is that, 
keeping the validity of (3.3), such an extra term yields the ergodicity. Note that 
the generator (3.8) is a special case of the one discussed in Chapter 9 of [17]. (See 
(9.25) combined with (7.12) there for an expression of the generator.) In particular, a 
unique solution to the martingale problem for C a ^ m defines an Ai (F)-valued Markov 
process, which henceforth we call the £ Q m -process. Intuitively, because of absence of 
the 'motion process', the law of this process is considered as continuum convolution 
of the continuous-state branching process with immigration (CBI-process) studied in 
[15]. (See (3.11) below.) In addition, Example 1.1 and Theorem 2.3 in [15] concern 
the one-dimensional version of the £ aim -process without the drift. The latter proved 
that the offspring distribution and the distribution associated with immigration of 
the approximating branching processes may have probability generating functions of 
the form s + c(l — s) a+1 and 1 — d(l — s) a , respectively. 

Proposition 3.3 Let m G M.(E)° . Then C a , m in (3.8) and A a . m in (3.1) together 
satisfy (3.3) with C = T(a + 2) and \t>(??) = ^(^(E) -1 ^) for any $ G Fq. Moreover, 
the C a . m -process has a unique stationary distribution Q a ,m with Laplace functional 

I Q a , m (dv)e- {l >> f) = e-^ 1+ ^\ feB + (E). (3.9) 

JM(E)° 

A random measure with law Q a , m may be called a Linnik random measure since it is 
an infinite-dimensional analogue of the random variable with law sometimes referred 
to as a (non-symmetric) Linnik distribution, whose Laplace transform appeared al- 
ready in (2.14). It is obtained by subordinating to an a-stable subordinator by a 
gamma process. (See e.g. Example 30.8 in [19].) Namely, letting {Y a (t) : t > 0} and 
{l(t) '■ t > 0} be independent Levy processes such that 

E[e- xy °®] = e- txa and E[e~ x ^} = e~ n °^ 1+x \ t, A > 0, 

we have for each c > 

E[ e - XY °M C »] = E[e~^ xa ] = e ~ clo ^ 1+xa) , A > 0. 

The first equality implies that 

P(F a ( 7 (c)) G ■) = / P( 7 (c) G dt)P(Y a (t) G •)• 
Jo 



13 



(3.9) clearly shows an analogous structure underlying, i.e., 

Qa,m(-) = / Qm(dT])Q aj r,(-), 

Jm(e)° 

where Q m is the law of the standard gamma process on (E,m). (See Definition 5 in 
[22]). It is also obvious from (3.9) that, as a t 1, Q a ,m converges to Q m . In addition, 
one can see that 

Hm£ Q , m *fa) = (77, — ) - (77, — ) + (m, — ) =: £^(77) 

for 'nice' functions ^, where f^f (0 = ^^(v + e ^r) ■ This is a special case of the 
generator of MBI-processes discussed in Section 3 of [21]. It has been proved there 
that Q m is a reversible stationary distribution of the process associated with C m . 

Proof of Proposition 3.3. As already remarked, if the term — a -1 (77, ^S) in (3.8) 
would vanish, (3.3) can be shown by essentially the same calculations as in the proof 
of Lemma 17 in [13]. (In fact, the change of variable z =: i](E)u/(l—u) in the integrals 
with respect to dz in (3.8) almost suffices for our purpose.) So, for the proof of (3.3), 
we only need to observe that (77, = for ^ of the form ^(77) = <^(r](EY l 7]) with 
$ G J-q. But this is readily done by giving a specific form of $. Indeed, for = 
(/!, /1) • • • (/i, f n ) the function \l/ takes the form \I> (77) = (77, /1) • • • (77, / n )(?7, l)~ n , from 
which it follows that 

After integrating with respect to r)(dr), the numerator on the right side vanishes. 

The argument regarding ergodicity is based on a well-known formula for Laplace 
functionals of transition functions. (See (9.18) in [17] for a much more general case 
than ours.) To write it down, we need only auxiliary functions called ^/-semigroup 
[15] because there is no 'motion process'. These functions form a one-parameter 
family {ip(t, -)}t>o of non-negative functions on [0, 00) and are determined by the 
equation 

^(t, A) = --^(t, \) 1+a - -4>(t, A), m A) = A (3.10) 
at a a 

with A > being arbitrary. An explicit expression is found in Example 3.1 of [17]: 



[l + (l-e-')A Q ] 1/a ' 



Let {?7 t : t > 0} be an £ Q)m -process, and for each 77 e M(E) denote by E v the 
expectation with respect to {r] t : t> 0} starting at 77. Then for any / G B + {E) and 
t > 

-<vt,f)] 



E v 



= exp 



■(77, V t f)- I (m,(V s fr)ds 
Jo 



(3.11) 
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where V t f(r) = ip(t, f(r)). As t — > oo the right side converges to 



exp 



roo 

- (m,(Vjr)dt = exp [-<m,log(l + /"))] 
Jo J 



since by (3.10) 



d_ 

di 



\og(l + (V t f(r)) a ) = -(V t f(r)r. 



This shows the ergodicity required and completes the proof. ■ 

Proposition 3.4 Suppose that m(E) > 1 and let Q a , m be as in Proposition 3.3. 
Then 

E Q a , m [ri(E)~ a ] = (T(a + l){m{E) - l))" 1 . 

Moreover, 

P«,M = r (« + - l)^ a ' m [v(E)- a ; viE)-^ e ■] (3.12) 

is a stationary distribution of the A a;m -proeess . 

Proof. The first assertion is shown by using t~ a = r(a) _1 J °° dvv a ~ 1 e~ vt (t > 0) 
and (3.9) with / = v. Indeed, these equalities together with Fubini's theorem yield 



r(a) 



r(a + i) 



dvv a ~ l exp [-m(E) log(l + v a )\ 



dzexp [-m(E) log(l + z) 



= r(a + l) -1 (m(£) -l)" 1 . 

As in the one-dimensional case, Theorem 9.17 in Chapter 4 of [7] reduces the 
proof of stationarity of (3.12) with respect to A a , m to showing that 



/ P a ^{dn)A a ^{n) = 

JMi(E) 



(3.13) 



for any $ of the form $(/i) = (/i, f\) • • • (/i, f n ) with /j G C(E') and n being a positive 
integer. Without any loss of generality we can assume that < fi(x) < 1 for any 
x G E and i = 1, . . . , n. Furthermore, we only have to consider the case where f\ = 
• • • — f n —'■ f because the coefficients of the monomial t± ■ ■ ■ t n in (//, tif\ + - ■ ■+t n f n ) n 
equals n\(/j,, f\) ■ ■ ■ (fi, f n ). Thus, we let $(/x) = (/i, f) n with < f(x) < 1 for any 
x <E E. Because of the basic relation (3.3) and (3.12) together, (3.13) can be rewritten 
as 

/ Q a , m (dri)C a , m *(ri) = 0, (3.14) 

JM(E)° 

where ^(rj) = (rj , /)" (rj , l)~ n . The main difficulty comes from the fact that \I/ does 
not belong to T . For each e > 0, introduce ^ e (v) := {"Hi /) n (( r ?> 1) +e) _ ™ and observe 
that \t> e G J 7 . Thanks to Proposition 3.3, we then have (3.14) with \& e in place of \I> 
provided that £ Qjm ^ e is bounded. Thus, the proof of (3.14) reduces to showing the 
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following two assertions: 

(i) For every e > 0, C^ m ^ € , C^ m ^ e and £^m*« are bounded functions on M(E). 

(ii) It holds that for each k G {1, 2, 3} 

"SfjL. = (3.15) 

Here, £ a , m = C^ m + C^ m + C^} m , and the operators £« m , £g) m and £gL correspond 
respectively to the first, second and last term on the right side of (3.8). 
First, we consider C^ m . Observe that 



Sv «77,l) + e)» ((r ? ,l) + e)«+ 1 

n(/(r)( T y,l)-( T y > /) + €/(r))(7y > /) B - 1 
((77, 1) + c)«+i 

from which it follows that 



(3.16) 



= —ne 



5r) 

n(( V ,f)( V ,l)-( V ,f)( V ,l) + e( V ,f))( V ,f) 
((7 7 ,l) + e)-+ 1 

^(77) 



n-l 



fa, 1) + e 



Hence £^ m ^ £ is a bounded function on AI(E') and £^^(77) ->■ = ££^(77) 
boundedly as e 4, 0. This proves that (i) and (ii) hold true for C^ m . 

In calculating £^^ e , (3.16) is useful since ^^(77 + z5 r ) = g^jj^ Indeed, 
by Fubini's theorem 

dz rT , „ , T , r°° dz r z , S^ f , , 

aw 



5(77 + u>5 r ) 

- / ^" Q ^T7 e ~T^(r), (3.17) 

a Jo 5(77 + w5 r ) K ' 



and combining with (3.16) yields 
/„ JI+J [*«(>! + z*0 - *«(•!)] 

Q! JO 



< - / °° ^.-^. n l/( r )^ + W(5 '-' !> - + /> + e/(r)| (77 + w8 r , f) n 1 



((77 + w5 r , 1) + e) n+1 



71 

< — I w~ a dw 



a Jo (77, 1) + w + e 



JO JO 



n 

= — I w 

a Jo 

= n r(a)r(1 " a) «»M) + e)-°. (3.18) 
a 
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This shows not only that is bounded but also 



\A*e(v)\<TiT{a) 



(to, 1) 



which is integrable with respect to Q a ,m as proved already. It can be seen also from 
(3.16) and (3.17) that C^ m ^ e converges pointwise to C^ m ^ as e | 0. By Lebesgue's 
dominated convergence theorem we have proved (3.15) for C^ m . 
The final task is to deal with C^ m ^ e . Similarly to (3.17) 



I e (r),r) 



^ dz 

y2+a 



J 

Jo 

r°° dz r 
Jo z 2+a Jo 



y e (r) + z5 r )-V e (r))-z 



5%, 
5rj 



dw 



— I 

1 + a Jo 



00 dw 



w 



l+a 



6V f 



5(rj + w5 r ) 



5^f 



By (3-16) ^y(r) - ^(r) equals 
((77, 1) + e) n+1 n (f(r) (77, 1) - (77, f) + ef(r)) Urj + w5 r , /) n_1 - (77, fY' 1 



+ 



((77, l)+w + e) n+1 ((r], 1) + e)^ 1 
((77, 1) + e) n+l - ((77, l)+w + e)" +1 l n (f(r)( V , 1) - (77, /) + e/(r)) (77, /) 



n-l 



((77, l)+w + e) n+1 ((r], 1) + e) n+1 
Moreover, we have bounds 

(77 + wd r , f) n ~ l - (77, f) n ~ l 



dv(n- l)f(r){r) + v5 r ,f) 



n-2 



< w(n — 1)((?7, 1) + w) 



n-2 



and 



((7 ? ,i)+ e r +i -((77,i)+ w + e ) 



n+l 



(n+1) dv((r), 1) + v + e) n 



< w(n + 1) ((77, l) + w + e) n . 



Consequently 

5(77 + w<5 r ) lr j 577 lrj 



n((77,l) + e) n+2 (n-l)((77,l)+7i;) 



n-2 



< W 



((77, 1) + w + e)" +1 ((?7, 1) + e) n+l 
{n + 1)((77, 1)+w + e) n n((v, 1) + e)(v, 



((77, l) + w + e) n+1 ((*7, 1) + e)" +1 
2n 2 



((77, 1) +w + e)(<77, 1) + e) 
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Therefore, analogous calculations to those in (3.18) lead to 

a + 1 



r(l-a) 7b 
< 2n 2 r(a)((r ? ,l) + e) 



I e (r),r)r)(dr) 



<?7, 1> + e" 



This makes it possible to argue as in the case of £^^ e to verify (i) and (ii) for C^ m . 
We complete the proof of Proposition 3.4. ■ 

Next, we show the coincidence of two distributions (3.4) (or (3.12)) and (3.7). Before 
going to the proof, it is worth noting that 



JMi(E) 



(3.19) 



where in general, for 6 > —a and m G A4(E), T>fy e ^ is the law of the two-parameter 
generalization of the Dirichlet random measure with parameter (a, 6) and parameter 
measure m defined by 

(See e.g. Section 5 of [22].) We will make use of the identity 

/ V^\diM)(jjL,l + f)- a =(m,(l + f) a )- 1 , feB + (E). (3.20) 
This is a special case of Theorem 4 in [22] and can be shown as follows. 

JMi(E) L v ' 

= T(a + l)E Qa ' m [(r],l + f)- a 

roo 

= a dvv a - 1 exp[-v a {m,(l + f) a )\ 

= wi+zn- 1 . 

Lemma 3.5 If m(E) > 1, then P a , m in (3.12) coincides with P a , m in (3.7). 
Proof. It suffices to show that for any / e B + {E) 



Mi(E) 



P a , m (dfi){fJl,l + f)- a 



Mi(E) 



P a , m (d»){»,l + f)- a =-- Hf)- 



In view of (3.12), calculations similar to the proof of (3.20) show that 
(r^ + lXm^)-!))- 1 ^/) = [fa, 1 + /)-"" 
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roo 

T(a)- 1 / dvv*- 1 exp \-{m, log(l + v a (l + /) Q ))1 
•/ o 

T^ + l)- 1 / dzexp[-(m,log(l + z(l + /) Q ))] 



— r / du(l — u) 2 exp 
1) Jo 



u 



-<m,log(l + - (1 + /D) 

1 — it 



r(« + i) 

i— T /' d«(l - M ) m ^- 2 exp [-(m, log(l + «((1 + /)" - 1)))] 
+ 1 Jo 



r(« + i) 
1 

T(a + 1) 



f 1 d«(i - ^r^- 2 / vjdfib, i + «((i + rr - i)))- m(E \ 

JO JMi(E) 



where the last equality follows from (3.6). Hence, by applying Fubini's theorem and 
(2.4) 



1(f) 



V m (d/j,) 



Bi :m (E)-l(du) 



i Mm JO (//, 1 + «((1 + /)« - i))) m ( s ) 

= / v m (d^,(i + fr)-\ 

JMi(E) 

On the other hand, combining (3.19) with (3.20), we get 

/(/) = / v^d^ii+frr 1 (3.2i) 

JMi(E) 

and therefore /(/) = /(/) as desired. ■ 

Remark. The 'semi-explicit' form (3.19) can be explicit if m is a probability measure. 
More precisely, we have P a>l/ = V au for any v G M.\(E). Indeed, observe that by 
(3.21) with m = v 

[ p a , u (dfx)(^i + fr a = I ^(^xMi+zrr 1 

JMi(E) JMi(E) 

= exp[-<i/,log{(l + /)"})] 
= exp[-(m/,log(l + /))] 

= / V av (d^,l + f)~ a , 

JMi(E) 

where (3.6) has been applied twice. (A one-dimensional version of the identity P a>l/ = 
V av is mentioned in Remark (ii) at the end of Section 2.) By (3.19) what we have 
just seen is rewritten as 

/ V u (dfi)V^\.)=V au (.), 

JMi(E) p 

which is a special case of 

/ vpfi'°\drivM>(-) = V^ e \-), (3 G [0, 1), 9 > -a(3. 

JMi(E) p 
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Here notice that, in case (3 = 0, T>^^ = T)g u by definition. This generalization can be 
proved analogously by virtue of the two-parameter generalization of (3.6) and (3.20). 
(See e.g. Theorem 4 in [22].) 

We can now prove our main result, Theorem 3.2. In the proof we write 9v {9 > 0, 
v G M.\(E)) for the parameter measure m. 

Proof of Theorem 3.2. Let v G M.\{E) be given. We first show that, for arbitrary 
9 > 0, P a ,8v is a stationary distribution of the ^^-process. For the same reason as 
in the proof of Proposition 3.4 (cf. (3.13)), it is sufficient to prove that 

/ PafiM^AafiMv) = (3.22) 

JMi(E) 

for <3> of the form $(/i) = (/i, f) n with / G C(E) and n being a positive integer. Since 
Proposition 3.4 and Lemma 3.5 together imply that (3.22) holds true for any 9 > 1, 
it is enough to show that the left side of (3.22) defines a real analytic function of 
9 > 0. We claim that 

A a ,eM») = ^ ± f) (1 - « W« + 1 h />-* - />") 

(3 (i -«)-(«)-«^*>^ ^ 

- ; E (") (1 - «)*- 2 (a + !)„_*</*, /*> fc, />""* 



T(n) fe f 2 V fc 

+7^W £ (^(l-a^a)^,/*)^/)-* (3.23) 



fc=i 

[a + l)„_i 



a+lfn 



(0 + n -!)</*,/>» 



The first equality is a special case of (3.2), and the second one can be shown with 
the help of Leibniz's formula 



(0i0 2 ) w (o) = E(^)#~ fc) (o)#(o) 



k=0 



for 0x(t) = (1 - t)' a and <p 2 (t) = (1 - t)" 6 with (a, 6) = (a + 1, -a - 1) or (a, 6) = 
(a, —a). In view of (3.23), it is clear that the proof reduces to verifying real analyticity 
of / P a ,eu(dfx)(fx, fi) ■■■ (/!, /„) in 9 for arbitrary /i, . . . , f n G C(-E). 

To this end, we shall exploit the following identity which is equivalent to (3.21): 

/ P a Md^(^l + f)- a =[ V ev (d^,(l + f) a )-\ (3.24) 

JMi(E) JMi(E) 

where / G B + {E) is arbitrary. Clearly this remains true for all bounded Borel 
functions / on E such that inf rG £ f(r) > —1. Therefore, for any t± : . . . , t n G R with 
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\ti\ + ••• + \t n \ being sufficiently small, (3.24) for / = — Y^i=\Ufi is valid, that is, 
J(*i, ••■,t n ) = J(t ir ■ ■ ,t n ), where 

(n \ ~ a 

l-(/^X>/i> (3.25) 
i=i J 



and 



(n \ a 

1-X>/* )"'■ ( 3 - 26 ) 
i=\ ) 

Noting that (1 — t)~ a = 1 + Y^k=i( a )kt k /k\ as long as \t\ is small enough, we see from 
(3.25) that the coefficient of the monomial t\ ■ ■ -t n in the expansion of I(t\, . . . ,t n ) 
is given by 

(«)n / P a! v (dll)(lI,f 1 )---(fJ,,f n ). (3.27) 

JMi(E) 

To find the corresponding coefficient for J(ti, . . . , t n ), define 



oo 



Kit) = 1 - (1 - t) a = a 2(1 - a),_i*7Z! 
and observe from (3.26) that J(ti, . . . ,t n ) equals 

Jm 1 {e) \ i=l J 



oo 

= 1 + 



k= l J M,{E) \ l=l J 
oo „ oofef/i _A / n \h • 



One can see that the coefficient of the monomial ti • • • t n in the expansion of J(ti, . . . , t n ) 
can be expressed as 



k=l T 67r(n,fc) ll J j=l ^ 1 Jjl ' i&lj 



(3.28) 



where 7r(n,k) is the set of partitions 7 of {l,...,n} into /c unordered nonempty 
subsets 71, • • • , 7fe. By Lemma 2.2 of [6] (or equivalently by Lemma 2.4 of [9]), each 
integral in the above sum is a real analytic function of 6 > 0. Hence, so is the integral 
in (3.27) and the stationarity of P a ,eu with respect to A a ,6v follows. 

It remains to prove the uniqueness of stationary distribution P of the A at e u -process 
for each 9 > 0. But this is an immediate consequence of (3.22) with P in place of 
P a ,6u and (3.23), which together determine uniquely / P(dfi)(fi, f) n and hence the 
nth moment measure 



M n (ciri • • • dr n ) := / P(d/j,)n(dr 1 ) ■ ■ ■ n(d 

JMi(E) 



(E) 
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for any n — 1,2, 



This completes the proof of Theorem 3.2. 



It is not clear whether we can derive from (3.28) an extension of the Ewens sam- 
pling formula in some explicit and informative form. (See Remarks after the proof 
of Lemma 2.2 in [6].) In view of (3.19), one might think that Pitman's sampling 
formula would be applicable. But it is not the case since V m (fx is discrete) = 1. The 
expression (3.12) might be rather useful for such a purpose. 

4 Irreversibility 

In this section we discuss reversibility of our processes. In contrast with the Fleming- 
Viot diffusion case, we guess that for any < a < 1 and non-degenerate m the 
A^m-process would be irreversible. Unfortunately, the following result does not give 
an affirmative answer in all cases. However, this does not suggest any possibility of 
the reversibility in the exceptional case, which is believed to be dealt with a different 
choice of test functions. 

Theorem 4.1 Let m e M.(E)° be given. Assume that either of the following two 
conditions holds. 

(i) The support of m has at least three distinct points. 

(ii) The support of m has exactly two points, say r\ and r2, and m({ri}) ^ m({r2}). 
Then the stationary distribution P a , m of the A a ^ m -process is not a reversible distri- 
bution of it. 

Proof. As in the proof of Theorem 3.2, we write 6u instead of m. Thus, 9 > and 
v E M.i(E). Recall that an equivalent condition to the reversibility of P a ,6v with 
respect to A a ,eu is the symmetry 

E [<M aA <t>'] = E WA a , e ^] , d>, e F , 

in which E[-j stands for the expectation with respect to P a ,eu- (See the proof of 
Theorem 2.3 in [6].) In the rest of the proof we suppress the suffix 'a, #z/ for simplicity. 
Let / G C(E) be given and define $ n (/f) = (/x, /)" for each positive integer n. We 
are going to calculate 

A := E [$ 2 -4$i] - E [$iA$ 2 ] • (4.1) 
For this purpose, observe from (3.23) that 

^i(^) = -^t(^/>-<^/», (4-2) 
a + 1 

a$m = (n,f) + ^-M m, f) + ^¥-(v,f 2 )-(e + i)(v,f) 2 (4.3) 

a + 1 a + 1 



22 



and 



3(a + l)(/x,/ 2 > </.,/> + (!- a) </.,/ 3 > 



e 



a + 1 
9 



3a(a + l)(v,f)(fi,f) 2 + 







a + 1 



3(l-a)a(z/,/ 2 )</i,/> 



. (1 - a) (2 - «)<!/, / 3 > - (a + 2)(0 + 2)</x, /)< 
a + 1 



(4.4) 



Combining (4.2) with the stationarity £'[^4$i] = 0, we get E[(fi, /)] = (u, /). There- 
fore, it is possible to deduce from (4.3) and -E^4$ 2 ] = 



+ i ) B[fc,/n = -?5L(„, / )» + (i + < 1 - a > 



a + 1 



a + 1 



f (^/ 2 >- 



Moreover, this equality between quadratic forms is enough to imply the one between 
symmetric bilinear forms: 

(9 + l)E[fa f) (/i, g)} = ^-(u, f) (u, g) + (l + ^—^-o) (u, fg), (4.5) 



a + 1 



a + 1 



where g G C(E) is also arbitrary. In the rest of the proof we assume that (u, f) = 0. 
This makes the calculations below considerably simple. By (4.5) 



(4.6) 



The equality -E[^4$ 3 ] = together with (4.4) implies that 

(a + 2){9 + 2)E[{/m, f} 3 } = 3(a + l)M lj2 + (1 - a) (l + ^fl) / 3 >- (4-7) 

These preliminaries help us calculate A in (4.1) as follows. By (4.3) and (4.4) 

9 



E 



a + 1 



E 



(/,,/) ((/,, f) -{9 + !)<//, /) 2 )) 



(a + 1) + a9 



E[{^ /> 3 ]- M 1>2 



a + 1 

and hence (4.7) yields 

(a + l)(a + 2)(0 + 2)A 
= [(a + 1) + a0] ' 

-(a + l)(a + 2)(0 + 2)M 1;2 



3(a + 1)M 1>2 + (1 - a) ( 1 + ^0 ) (z/, / 3 ) 



a + 1 



= (a + l)(a - 1)(2(9 + l)M lj2 + [(a + 1) + a9\ (1 - a) 1 + 



2 - a 
a + 1 
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Plugging (4.6) into this expression, we obtain 

1 — a 

(a + l)(a + 2)(9 + 2)A = (q + + ^ fo f), 

where 

U(a,6) = -(a + l)(20 + l)[(a + l) + (l-a)0] 

+ [(a + 1) + a9\ (9 + 1) [(a + 1) + (2 - a)9\ 
= a9 2 [(a + 4) + (2 - a)9] =: V(a,9). 

(The second equality between quadratic functions of a is verified by checking that 
U(-l, 9) = -39 2 (9 + 1) = V{-\, 6), U(0, 9) = = V(0, 6) and 17(1, 0) = 9 2 (9 + 5) = 
V(1,0).) Consequently, whenever (v, f) = 0, we have 

a(l - a)9 2 [(a + 4) + (2 - a)9] 3 



(« + l) 2 (a + 2)(^ + l)(^ + 2) 



Thus, all that remains is to construct an / G C(E) such that (u, f) = and 
(u, f 3 ) > 0. Because of the assumption, we can choose a closed subset E of E such 
that < v(Eq) < 1/2. Indeed, in the case (ii) this is trivial while in the case (i) there 
exist disjoint closed subsets Ei,E 2 and E 3 of E such that v(Ei)v(E 2 )i>{Ez) > and 
so < v(Ei) < 1/2 for some % G {1, 2, 3}. Letting g denote the indicator function of 
E Q , we observe that 

(u,(g-(u,g)) 3 ) = (v : g 3 )-3(v,g 2 )(v : g) + 3(v,g)(v : g) 2 -(v : g) 3 
= u(Eo) - 3u(E ) 2 + 2u(E ) 3 
= u(E )(l-u(E ))(l-2u(E ))>0. 

Finally, the required / exists since g can be approximated boundedly and pointwise 
by a sequence of functions in C(E). The proof of the theorem is complete. ■ 

It is worth noting that the exceptional case of Theorem 4.1 corresponds to a subclass 
of the one- dimensional case discussed in Section 1, more specifically, the process gen- 
erated by (1.3) with c\ = c 2 . There is no reason why this class should be so special 
with respect to the reversibility, and it seems that such a 'spatial symmetry' makes 
it more subtle to see the asymmetry in time. The actual difficulty in showing the 
irreversibility for these processes along similar lines to the above proof is that expres- 
sions of E[& ril A& n2 ] with ni + n 2 > 4 as functions of a and 9 are too complicated to 
handle. 
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